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Abstract 

This work is an extended version of a paper published last year in Physical Review Letters 
[S. de Man et al, Phys. Rev. Lett. 103, 040402 (2009)], where we presented measurements of the 
Casimir force between a gold coated sphere and a plate coated with either gold or an indium-tin- 
oxide (ITO) layer. The experiment, which was performed in air, showed that ITO is sufficiently 
conducting to prevent charge accumulation, but still transparent enough to halve the Casimir 
attraction when compared to gold. Here, we report all the experimental details that, due to the 
limited space available, were omitted in the previous article. We discuss the performance of our 
setup in terms of stability of the calibration procedure and reproducibility of the Casimir force 
measurement. We also introduce and demonstrate a new technique to obtain the spring constant 
of our force sensor. Furthermore, we present a thorough description of the experimental method, a 
comprehensive explanation of data elaboration and error analysis, and a complete characterization 
of the dielectric function and of the surface roughness of the samples used in the actual experiment. 
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INTRODUCTION 



It is well known that the Casimir effect [T] strongly depends on the dielectric function of 
the interacting surfaces [21 E]- Transparent dielectrics, for example, attract less than highly 
reflective metals. Dielectric materials, however, tend to accumulate isolated charges. Those 
charges give rise to an electrostatic force that easily overcomes the Casimir interaction. 

In a recent paper [1], we have presented measurements of the Casimir force between 
a gold coated sphere and a plate coated with either gold or an Indium-Tin-Oxide (ITO, 
In203:Sn) layer. The experiment, which was performed in air, showed that ITO is sufficiently 
conducting to prevent charge accumulation, but still transparent enough to halve the Casimir 
attraction when compared to gold. 

The experiment was carried out by means of a quite complicated novel technique that, 
due to the hmited space available, was not thoroughly explained in our previous work. We 
believe it is important to extend that work and provide the community with all the details 
of the experimental technique and data analysis, which is the purpose of this paper. 

This paper is organized as follows. First, we describe the experimental setup and discuss 
general issues one has to tackle to perform Casimir force measurements. Then we discuss 
the experimental technique we developed to simultaneously calibrate the setup and measure 
the Casimir force gradient, and derive in detail the specific forms of all our calibration and 
measurement signals. Second, we illustrate a new method to determine the spring constant 
of our force sensor. Third, we present experimental results on the general performance 
of our setup, namely the stability of the calibration procedure, the reproducibility of the 
force gradient measurements, and the spring constant determination. Fourth, we present 
the Casimir force measurements for the gold-gold and gold-ITO interactions, and show 
measurements of the dielectric functions and surface topographies of our surfaces. Finally 
we compare the hydrodynamic forces for the two sets of experiments. 

EXPERIMENTAL SETUP 
Description 

Our experimental setup is designed to precisely measure surface forces between a 100 fim 
radius sphere and planar samples at ambient pressure. The sphere is attached to a mi- 
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cromachined cantilever (spring constant roughly 1 N/m) whose deflection in response to 
external forces can be measured with pm sensitivity by a commercial Atomic Force Micro- 
scope (AFM) detection head (Veeco Multimode); the detection system is formed by a laser 
beam that reflects from the free end of the cantilever and hits a position sensitive photode- 
tector (see Fig. [T]). The sphere-cantilever assembly is coated with a Ti adhesion layer and a 
100 nm Au film. The planar sample is mounted on a two-stage mechanical translator formed 
by a stick-slip piezoelectric motor (Attocube) and a feedback controlled piezoelectric trans- 
ducer (Physik Instrumente) to vary the separation between the sphere and plate surfaces. 
The stick-slip motor is used for coarse approach (travel range 6mm), while the feedback con- 
trolled transducer executes the fine distance scanning (range 12 fim, closed loop resolution 
50 pm). Both the detection head and the two-stage mechanical translator are anchored to a 
10 cm^ Al block that is actively temperature stabilized at 300 K to reduce mechanical drift 
from differential thermal expansion of the components. The Al block is screwed onto an 
active ant i- vibration table (Halcyonics) , which is placed inside an anechoic chamber. This 
chamber lies onto a heavy marble optical table that is located in a temperature controlled 
laboratory. 

Three crucial issues 

In a Casimir force measurement, there are three crucial issues that have to be dealt with. 

First, even if one would electrically connect both interacting surfaces, there exists an 
electrostatic potential difference Vq due to the different work functions of the surfaces. Since 
work functions of surfaces depend on quite a number of parameters, like crystal growth 
orientation and adsorbates, typically there even exists a potential difference between surfaces 
made out of the same material. This electrostatic potential difference gives rise to a force 
that is generally stronger than the Casimir force. To avoid this problem, most Casimir force 
setups rely on a counterbias circuit that is used to apply — Vq to the surfaces in order to 
have no residual electrostatic force. 

Second, even in setups where the distance between the sphere and the plate is varied with 
a feedback controlled piezoelectric transducer, one has only knowledge of the relative position 
changes and not of the absolute separation between the surfaces. It is thus mandatory to 
find the initial separation do with a calibration procedure. Because the distance dependence 



3 



of the electrostatic force between a sphere and a plate is known exactly, most modern setups 
use this force to extract cIq. 

Third, the instrument has to be calibrated with a known force. Again, one can use the 
electrostatic interaction to calibrate photodetector voltage versus force. We have developed 
a measurement scheme that solves all three issues at the same time. 

Force modulation measurements 

We present a measurement technique that makes use of simultaneous detection of both 
calibration signals (based on the electrostatic force) and the Casimir force. The motivation 
for this approach is the benefit of absolute certainty that the calibration parameters always 
correspond to the measured forces because they are acquired simultaneously; it is thus 
impossible to have inconsistent calibration and force data due to time-related drifts or other 
events. In order to achieve this goal, we have separated the calibration and Casimir signals in 
frequency space: the signals are modulated at distinct frequencies that can be de-modulated 
individually with lock-in amplifiers. 

Modulating an electrostatic interaction is extremely easy: one just has to apply a time- 
dependent potential difference to the sphere and the plate {Vac in Fig- [l]b). We thus apply 
an oscillating voltage Voc + Vac cos (wit) between the sphere and the plate, where V^c is 
used to compensate for the contact potential difference Vq [5] • Unfortunately, modulating the 
Casimir force is a lot more challenging as its strength depends only on geometry and dielectric 
properties of the surfaces. On the other hand, of course, the strong distance dependence 
of the Casimir force can be used to modulate its strength considerably. Therefore, we 
chose to add a small modulation of the form Ad cos (^2^) to the piezoelectric transducer 
displacement dp^, as previously introduced in [S]. When the sphere and plate surfaces 
are separated by a distance d, we have the following three forces acting on the sphere: 
F{V, d, U2, Ad) = -Fe(V, d) + Fc{d) + Fn^d, U2, Ad) where Fe^V, d) is the electrostatic force 
for externally applied potential difference V, Fc{d) is the Casimir force, and FH{d,U2, Ad) 
is the hydrodynamic force due to the moving air caused by the oscillatory motion of the 
plate. These forces induce a bending of the cantilever F/k according to Hooke's law, where 
k is the spring constant of the cantilever. The output of the optical lever read-out S is 
then changed by AS* = 'jF/k, where the sensitivity of the read-out is characterized by the 
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calibration factor 7. We will now develop the full form of this signal AS. 

Following elementary electrostatic arguments, one can show that the electrostatic force 
between a plane and a sphere of radius R is given by 

where Eq is the permittivity of vacuum, V is the externally applied voltage, Vq is the contact 
potential difference between the two surfaces, and d <^ R {i. e. within the proximity force 
approximation (PFA) P]). To evaluate the total signal AS", we substitute V = Vnc + 
Vac cos (uit) and incorporate the distance modulation Adcos{u2t). We then approximate 
the photodetector signal with a first order Taylor expansion for small excursion Ad cos {u)2t) 
around d = do — dpz (see Fig. [T] c): 



AS{t) ~ ^0 + S^, cos{uJit) + cos(2a;it) + S^^ cos{uj2t) + 5*^^ sm{u2t) + 5rem(t), (2) 
where 

^° k (do - d,z) + fc^^ ~ ^^^^ ' 

27eo7ri? (Vp + Vpc) Vac , . 

k{do-d,.) ' 



, ^ 7govri? [(K) + VpcY + ylciA 7 dFc 

k{do-d,zf k dd 



Ad, (6) 

do— dpz 



SZ = lFH{do-dp,,U2,Ad), (7) 

-^c {do — dpz) is the Casimir force at separation do — dpz, Stemit) contains the cross terms at 
frequencies like Ui ± UJ2 and 2uJi ± 002 and the gradient of the hydrodynamic force at 2uj2, 
and we have neglected the effect of the cantilever deflection on the distance between the 
surfaces. Since the remaining terms in Srem are located at different frequencies than our 
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measurement signals, they do not interfere with the lock-in measurements of S^jj^, S2UJ1, S^^ 
and • ^rcm will thus be neglected in the rest of the paper. 

Eq. [2] is only valid if the force sensor can follow the modulations of the force without 
picking up phase delays. It is thus convenient to operate in the quasi-static regime, which also 
ensures that the amplitude response of our cantilever at the various measurement frequencies 
does not vary. For these reasons, we set u)i/2tt = 72.2 Hz and u)2/2tt = 119 Hz, which are 
both much lower than the resonance frequency of the force sensor (1.9 kHz, quality factor 
~ 75 in air). Furthermore, we have not included the elastic component of the hydrodynamic 
interaction in S^^. According to [7], the compression effect is small as long as a sphere = 
^^^^ < 1, where rj is the viscosity of air and p is the air pressure. In our experiment, 
o'sphere < 10"'^, SO the clastic component can be neglected and we only have to consider the 
dissipative part of FH{d,uj2). Since a dissipative effect depends on velocity v = dd/dt = 
U2Ad sin{uj2t) , it will manifest itself as a cantilever oscillation at 002 with a corresponding 
detector signal S^^ that is 90 degrees rotated with respect to S^^. 

Electrostatic calibration 

The first task of the electrostatic calibration procedure is to compensate for the presence 
of the contact potential difference Vq between the two interacting surfaces. S ince Sijj. is 
proportional to Vq + Vdc (see Eq. |4]), we can create a negative feedback loop in which a 
lock-in amplifier at coi generates Vdc in such a way that S^j-^ vanishes, i.e. Vdc = —Vq 
O E]. The stability of this feedback loop is guaranteed by a single large time constant. 
In the current experiment, the systematic error in the compensation voltage is negligible 
(|Vo + Vdc\ < 50 yuV), and the statistical error is ~ 1 mV. This feedback scheme is similar 
to Kelvin probe force microscopy [9j, and allows one to measure Vq at all sphere-plane 
separations. Even more, the automatic compensation of Vq leads to the zeroing of the 
(Vo + Vdc) terms in Eqs. |3]and[6| greatly simplifying the measurement scheme. 

The periodic component of AS* at 2uJi, 821^-^^, measured with a second lock-in amplifier 
(the calibration lock-in in Fig. [l]b), is used to calibrate the force sensitivity and to find the 
initial separation between the surfaces do. We define 



According to Eq. [sj a can be experimentally obtained from a = 2\S2loi \ /V^^. ^^is 
way, we have essentially performed an AC measurement of the curvature of the electrostatic 
parabola, instead of using multiple DC measurements with different applied voltages [TU] . 
We measure a as a function of dpz by varying the extension of the capacitive feedback 
controlled piezoelectric transducer (see Fig. [T]:) in discrete steps. We then fit Eq. |8] to 
these a data. This procedure allows us to calibrate the separation at the start of the 
measurement and the force sensitivity k, = 'jeQiiR/k for each measurement run. We then 
use the estimate of do to adjust the initial value of dp^ of the next measurement run in order 
to have all runs start at the same separation. To avoid large electrostatic forces at small 
separations, we reduce Vac as the surfaces approach such that S2UJ1 stays nearly constant at 
a value that corresponds to a root-mean-square electrostatic force of ~ 50 pN [5] . 



Casimir force measurement 

We use a third lock-in amplifier (the measurement lock-in in Fig. [l]b), locked at UJ2, 
to measure the Casimir force. The phase of this lock-in amplifier is aligned to the actual 
motion of the plate by examining the signal from a dedicated fiber optic interferometer 
(not shown in Fig. [T]). The same interferometer is used to calibrate the amplitude Ad of 
the separation modulation. We see from Eq. [6] that the in-phase component S^^ contains 
both an electrostatic contribution and the gradient of the Casimir force Fc at the current 
separation. Since Vq + Vdc = by the Vq feedback circuit, Eq. [6] simplifies to 



k{do-dp,y 2 k dd 

Combining Eqs. |5]and[9| one obtains 
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Ad. (9) 



^ '5*2^^1 _ 7 dFc 
'^^ do — dpz k dd 



Ad. (10) 

da— dpz 

Since the absolute separations do— dpz and S'2a;i are known from the simultaneous electrostatic 



calibration (and Ad is calibrated too), one can calculate the value of the first term of Eq. 10 
Using the force sensitivity n = '-feoTcR/k obtained from the calibration, we can finally get 
the Casimir force gradient: 



It is interesting to note that we obtain the Casimir force gradient divided by the sphere 
radius R, because we have cahbrated the instrument with the electrostatic force which 
scales linearly in R (see Eq. [T|. However, within the PEA, the gradient of the force between 
a sphere and a plate relates directly to the pressure between two parallel plates Ppp as long 
as d R: 

i|f=2.P„M, (12) 

where Ppp{d) can be calculated with the Lifshitz theory [2J and depends only on the dielectric 
properties of the interacting surfaces. Therefore, we can directly compare our 1/R dFc/dd 
data to theory, without any need to know the precise radius of the sphere. 

Eurthermore, by using a quadrature lock-in amplifier at uj2, we can obtain S^^ together 
with 5*^2 . We can thus measure the hydrodynamic interaction between the sphere and the 
plate simultaneously with, but independently from, the Casimir force gradient. 



Determination of deflection sensitivity and cantilever spring constant 

So far, we have neglected the bending of the cantilever in the assessment of the distance 
between the sphere and plate surfaces. This is valid as long as the forces are relatively 
weak and the spring constant of the cantilever is relatively high. Of course, the nominal 
spring constant of the cantilever is supphed by the manufacturer, but the addition of a glued 
sphere and metal coating influence the stiffness. Therefore, we have developed a technique 
to measure the spring constant with the electrostatic force. Eurthermore, this method also 
allows us to extract the deflection sensitivity 7 of the optical lever readout; we can then 
convert photodetector signal AS* into cantilever deflection F/k. This technique might be 
useful for AEM force measurements in general. 

To obtain the cantilever spring constant and the deflection sensitivity, we apply a rel- 
atively large Vac between the sphere and the plate. Exactly like described above in the 
electrostatic calibration section, we keep the electrostatic force at 2uJi constant, but now at 
roughly 2 nN RMS instead of 50 pN RMS, by reducing Vac while increasing the piezoelec- 
tric transducer extension dpz in discrete steps. This strong force will reduce the sphere plate 
distance, and we therefore have to solve the following implicit equation for the electrostatic 
force 
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where, since we have already dealt with the contact potential difference Vq with the feedback 
circuit, V just refers to the AC component of the applied voltage. For the sake of simplicity, 
we have omitted the piezo modulation at U2 from this derivation, as it does not affect the 
results. Eq. [13] has two solutions for Fe, and the physically correct one reads 

1 



k{do - dp,) - sJk'^{dQ - dp,y 



AkstMiRV^ 



(14) 



If we Taylor expand this expression for small cantilever deflection (which means small force 
and small applied voltage V), and use AS = •yF/k, we obtain 

~ ~k{do-dp,) - k\d,-dp,)^ + ^^^^ 
Substituting V = Vac cos(a;it) and neglecting the higher order terms yields a detector signal 

AS{t) ~ ^0 + ^2..i cos(2wit) + S^^, cos(4wit), (16) 

where Sq is the DC component and the amplitudes of the two AC components are given by 



2k{do-dp,) 2k^{do-dp,j 



\3 

and 



8k^{do - d, 



'pz ) 



S2UJ1 is already measured by our electrostatic calibration lock-in amplifier, and we simply 
add another lock-in amplifier locked at iui to detect 

The second term in Eq. [17] is much smaller than the first term and can be neglected. We 
then find that 

^ (do-dpz), (19) 




eoTiR 



which means that we can obtain k/R hj fitting Eq. 19 to data of Vac^/S2uii/S4uji as a 
function of relative piezo displacement dp,. Apart from the resulting knowledge on the 
cantilever spring constant (the sphere radius is roughly known), we also obtain the deflection 
sensitivity 7 by combining the value of k/R with the one of k = '-feoiiR/k determined by the 
analysis of the simultaneously acquired a data (as described in the electrostatic calibration 
procedure). 



RESULTS AND DISCUSSION 



We have divided our experimental results into two parts. In the first part, we describe 
the precision and stability of the electrostatic calibration procedure and comment on the 
reproducibility of the Casimir force gradient detection in a set of 580 measurement runs 
between two gold coated surfaces. Also, we present a single dataset obtained with a large 
electrostatic force between the sphere and the plate that allows us to check the validity 



of Eq. |6] (and Eq. 11 as well), and to obtain the spring constant of our cantilever and the 
deflection sensitivity of the optical lever readout. In the second half, we combine the Casimir 
force measurements between the two gold surfaces with measurements between a gold surface 
and a surface coated with ITO (In203:Sn), as presented in adding details that, for the 
sake of brevity, were previously omitted. Furthermore, we obtain the hydrodynamic forces 
for both measurement sets and compare the results. 



General performance 

We will now analyze the 580 measurement runs between two gold surfaces obtained during 
nearly 72 hours of continuous data acquisition. In this experiment, the separation between 
the surfaces is varied in discrete steps with the feedback controlled piezoelectric transducer, 
and a typical measurement run consists of ~ 50 dpz set points in the measurement range 
50 < (i < 1100 nm. The lock-in measurements are obtained with 24 dB roll-off low-pass filter 
settings with 1 s RC time. The waiting time for every value of dpz is 8 s, and a complete 
run takes roughly 7 minutes. The 82^}^ set point corresponds to a cantilever movement 
of approximately 50 pm RMS. The distance modulation is set to Ac? = 3.85 ± 0.08 nm, 
and the in-phase and out-of-phase cantilever responses at uj2 are < 80 pm RMS during the 
entire experiment. All force measurements are performed in air at atmospheric pressure, 
temperature 300 K, and 29% relative humidity. 

Concerning the electrostatic calibration, we have to fit our a data with Eq. [8]to obtain the 
initial separation and the force sensitivity k. Due to the fact that we hold 82^01 constant 
by reducing Vaci the relative statistical error in a is constant (see reference [3]) and was 
measured to be ~ 0.7%. We have verified that a follows Eq. [sj as suggested in jTT] and 
discussed in j5j . In Fig. |2] and Fig. [3| we present the fitted values for do and n and analyze 
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their stability in time. Fig. [2] a shows all the values of do for the 580 runs with error bars as 
propagated from the error on a. The grey line represents the smooth thermal expansion of 
the setup, and is estimated by smoothing the data with a 100 point moving window second 
order Savitsky-Golay filter. The total mechanical drift of our setup is 52 nm in 72 hours, 
which is less than 1 nm per hour and less than 0.1 nm per measurement run. Clearly, we 
can neglect the mechanical drift in our assessment of separation between the surfaces in a 
single run. In Fig. [2] b, we plot a histogram of the difference between the do data and the 
smoothed grey line of Fig. |2] a. These differences are clearly normally distributed with a 
standard deviation of 0.5 nm. Therefore, in these experimental runs, we could determine 
the separation between the sphere and plate surfaces with 0.5 nm precision. This estimate 
of the precision in the measurement of do is insensitive to the precise form and size of the 
smoothing window. 

Fig. [3] a shows all the values of the force sensitivity k that we obtained from the fit to 
our electrostatic calibration data. The error bars are calculated by propagating the errors 
on a. The grey line is a smoothed trend line that represents slow variations in n over time, 
obtained by smoothing the data with a 200 point moving window second order Savitsky- 
Golay filter. There is clearly no long-term drift in the force sensitivity, which shows that our 
setup is very stable. In Fig. [3]b, we plot a histogram of the relative deviations between our 
K data and the smooth grey line of Fig. [3] a. These deviations are normally distributed with 
a standard deviation of 0.2%, and are insensitive to the specifics of the smoothing window. 
We have thus determined the force sensitivity of the setup for every single measurement run 
with a precision of 0.2%. 

In Fig. |4] we present measurements of the total force gradient 

IdF BoTTSi^ 



Rdd K Ad 

as a function of the non-modulated separation d = do — dpz (see Fig. [T]c). This force gradient 



should, according to Eq. 11, obey 



Rdd R dd R dd 
with 

R dd K, do — dpz 



IdF IdFc^ldFE 



(22) 
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In Fig. |4] a, the data points represent the —eQ7rS^^/{K,Ad) data points and the sohd hne 
shows the electrostatic force gradient as obtained with Eq. 22 from the 6*2^^1 values of the 
simultaneous electrostatic calibration procedure. For clarity, we have shown only 150 mea- 
surement runs out of the total 580. It is clear that in the distance range that our setup is 
sensitive for the Casimir force, the electrostatic force gradient caused by the simultaneous 
electrostatic calibration is small compared to the Casimir force gradient. To assess the sta- 
bility of our force gradient measurement, we have plotted all 580 1/R dF/dd data points 
gathered around 95 nm from our 580 measurement runs in Fig. |4]b. Our data do not show 
any drift in time, which means that the setup is stable. The grey line represents the average 
of the data points. In Fig. |4]c, we plot a histogram of all the relative deviations of the data 
with respect to the average. These deviations are normally distributed with a standard de- 
viation of 3.5%, which corresponds to a standard deviation in the measurement of the force 
gradient of 1.85 N/m^. This value represents an overestimate of the noise though, because 
the data are obtained at slightly different separations; the exact position of a data-point 
depends on the estimate of do coming from the previous measurement run. Since the error 
in the determination of do is 0.5 nm (see Fig. [2]b), the data are horizontally scattered with 
a standard deviation of 0.5 nm. For d ~ 95 nm, the local slope of the data in Fig. |4] a is 
approximately 1.3 Nm~^/nm, which translates this scatter in d into a force gradient scatter 
of 0.65 N/m^. Therefore, the actual precision in a single force gradient data point around 
95 nm is 1.75 N/m^, if we assume that both the force gradient noise and the scatter in 
separation are uncorrected. 

From the electrostatic calibration results, we could have also estimated the noise in the 
force gradient measurement. In fact, the noise in S2UJ1 is 30 /iV RMS with a 1 s RC time. 
The force gradient signal at U2 is located at a comparable frequency, therefore the noise will 
be quite the same. If we substitute our measured values of k and Ad into Eq. 20, we see 
that we would have expected the noise in 1/i? dF/dd to be 1.62 N/m^. But we have not 
taken into account yet the 0.5 nm error in the separation that arises from the estimate of 
dQ. For (i ~ 95 nm, this results in an additional statistical error of 0.65 N/m^ in the force 
gradient at this distance. The combined error, assuming the force gradient and distance 
errors are uncorrelated, is then 1.75 N/m^ for (i ~ 95 nm, which agrees perfectly with the 
data of Fig. |4) 

Since our Casimir force gradient measurement consists of measuring the total force gradi- 
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ent and subtracting the electrostatic force gradient (see Eq. 11 ), it is interesting to investigate 
the accuracy in the assessment of 1/i? dFE/dd. For that, we have gathered a new dataset 
with a relatively strong electrostatic interaction (high Vac) between the sphere and the 
plate. When we combine the high Vac total force gradients with measurements obtained 
with low Vac, we can get 



1 dF 
R~dd 



1 dF 
R~dd 



1 dFp 



1 dFf 



V^c> ^ 



(23) 

Vac< 



Vac> 

because the Casimir force gradient is equal in both cases and drops out. Even more, any other 
systematic effects present in the force gradient measurement that do not depend on Vac, like, 
for example, laser light that reflects from the planar sample and hits the photodetector, are 
also cancelled in this way. The right-hand-side of Eq. [23] can be calculated with Eq. |22| and 
we can thus assess the validity of the latter and, consequently, of EqjlTj In Fig. [5| we have 
plotted the difference in total force gradients (obtained with Eq.|20]) as a function of distance. 



The solid line represents the difference in calculated electrostatic force gradients (Eq. 22), 
determined with the corresponding sets of 6*2^^1 data. Although the agreement between the 
two electrostatic force gradients is good (there are no adjustable parameters), there exists 
a slight discrepancy between the two curves. The measured total force gradient difference 
is systematically about 3% higher than the values calculated from S2uj-^- If this discrepancy 
means that there is a small error in the determination of the electrostatic force gradient, 
then the measurement of the Casimir force gradient is almost unaffected. For example, for 
all d < 120 nm the electrostatic force gradient is always < 25% of the total force gradient, 
which results in an error of < 1% in the measurement of the Casimir force gradient. If, on 
the other hand, the mismatch is caused by the uncertainty in the determination of Ad with 
the dedicated fiber optic interferometer, then our Casimir force gradients are affected by a 
3% systematic error. Nevertheless, this systematic error will not hamper the comparison 
between force gradient data obtained with different samples, as we always use the same Ad. 

To measure the deflection sensitivity of the readout and the spring constant of the can- 
tilever, we follow the procedure outlined above. In essence, we apply a big potential difference 
between the sphere and the plate, record both the cantilever deflection signal at 2u)i and Aui 



as a function of relative piezo-electric transducer displacement dpz, and fit Eq. 19 to those 



data. In Fig. joj we plot Vac\/ S^uj^/ S^^^^ as a function of dpz for such a single dataset. The 



straight line represents the best fit with Eq. 19 (reduced x = 0.25). The error bars were de- 
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termined by measuring the absolute error in S2UJ1 and assuming that the error in 84^11 is equal 
and independent from the error in S2uji- This is not entirely correct, because some sources of 
error, like for example fluctuations in d, will lead to correlated variations in S2u]i and S^^^j^. 
We have thus overestimated the error in Vac \/ / S^i^^ , which leads to a reduced <\. 
Anyhow, the slope of the data allows us to extract k/ R= (11.12±0.06) 10^ N/m^ (the uncer- 
tainty is obtained by setting reduced = 1). When we combine this value oi k/R with the 
simultaneously determined k = 191.3 ± 0.2 nm/V, we find that 7 = (7.64 ± 0.04) 10'' V/m. 
With this value of 7, we can now establish that the S2U11 set point we used for this dataset 
corresponds to a cantilever motion of 2 nm RMS at 2ui. It is interesting to observe that 
this 2 nm modulation of the separation d at 2uJi gives rise to a measurable signal at 4ui 
even at 1 /im distance. Furthermore, if we use the approximately known sphere radius of 
100 fim, we obtain the spring constant of our cantilever k = 1.1 N/m. As the nominal spring 
constant before sphere attachment and gold coating was 0.9 N/m, the value we find with 
this electrostatic method is very reasonable. 

With the deflection sensitivity calibrated, we can now assess the total cantilever bending 
and the precision in the measurements of the cantilever deflection. In the measurement runs 
presented in Figs. |2| |3] and |4], we used an 82^1 set point of 4 mV RMS, which corresponds to 
a cantilever motion of 52 pm RMS. Therefore, the static bending of the cantilever due to the 
electrostatic calibration procedure is 74 pm (see Eqs. [sjandjs]). Anyway, this static bending 
is constant during the measurement run and it is thus automatically taken into account in 
the estimate of d^. The cantilever oscillations at U2 caused by the total force gradient and 
the hydrodynamic interaction are < 80 pm RMS for these measurement runs, which means 
that the corresponding static bending is < 113 pm. It is thus evident that we can safely 
neglect the static bending of the cantilever in our data analysis. Since the noise in >S'2a;i is 
30 /iV RMS, the precision in the detection of the cantilever deflection is 400 fm RMS with 
our 1 s RC time (24 dB low-pass filter). This means that our setup has an RMS sensitivity 
of 1 pm/yHz at 2uJi/2n = UAA Hz. 

Halving the Casimir force 

We now present a comparison between two experiments performed with the same gold 
coated 100 fim radius sphere and two different plates. The first experiment is conducted 
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with a polished sapphire substrate coated with a gold film similar to the one deposited on 
the sphere. The general performance of our setup was discussed above by analyzing this first 
experiment. The second experiment consists of 580 measurement runs in which the plate is 
replaced by a float glass substrate with a sputtered ITO thin film on top (PGO CECOIOS, 
typically 8.5 Q/O, or, equivalently, p = 1.6 10~^ r2cm). After purchase, this sample has 
been exposed to air for more than two years before our measurements were performed. 

Fig. [T] shows the Casimir force gradient between the two pairs of surfaces (Au-Au in green 
triangles, Au-ITO in red squares) [1]. In Fig. [7]a, we plot the force gradients as a function 
of separation on a double logarithmic scale for randomly chosen subsets of the data (150 
out of 580 for both cases). Both datasets are obtained with the exact same settings for the 
electrostatic calibration and the force gradient measurement, and the Casimir force gradient 
is obtained from Eq. [TT] The black lines indicate the theoretical force gradient, as will be 
explained below. Figs. [7] b and c present two histograms of all 580 Casimir force gradient 
measurements for both Au-Au and Au-ITO at separations d = 120 nm and d = 80 nm, 
respectively. It is clear that the interaction strength with the ITO sample is considerably 
reduced with respect to the gold plate. 

Note that our estimate of do, and thus d, relies on the simple form of Eq. [T]and is only 
valid foT d R (PFA). This assumption is not entirely correct in the probed separation 
range |5j and results in a systematic error in rfo of about 1.4 nm. Still, the corresponding 
underestimate of the separation is equal for both the measurements with Au and ITO, and 
can thus be neglected in the comparison of the two experiments. 

Concerning the compensation voltage, we observed that Vq varies approximately 1 mV and 
3 mV over the complete measurement range in the Au-Au and Au-ITO cases, respectively. 
These slight variations of Vq do not compromise the measurement of the Casimir force at the 
current level of sensitivity. The value of Vq drifts in time from —106 to —103 mV for Au-Au 
and from 72 to 50 mV for Au-ITO at d = 100 nm. It is also important to note that, during 
the whole duration of the experiment, we never observed any problem with electrostatic 
charging of the Au or ITO layers, which would have most likely resulted in erratic behavior 
of a and/or Vq. 

In order to compare the obtained Casimir force gradients with theoretical predictions, we 
have investigated the dielectric properties of our surfaces. In Fig. |8j we show the reflection 
and transmission spectra of the two plates, measured from the thin film side, in the frequency 
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range from 0.5 to 6.5 eV. The black continuous lines represent the reflection and transmission 
spectra calculated from the literature. For Au, we used the values reported in [12]. The 
imaginary part of the dielectric function of ITO is constructed from a sum of Drude and 
Tauc-Lorentz models with the parameters from [ISj. The real part of the dielectric function 
is calculated with direct Kramers-Kronig integration. The thickness of the ITO thin fllm is 
fltted by examining the interference fringes in the reflection and transmission spectra (taking 
into account the refractive index of the material) and turned out to be 190 nm, which is close 
to the typical thickness reported by the manufacturer (180 nm). The agreement between the 
spectroscopic measurements and these literature values is reasonable in the probed energy 
range. We want to stress that for Au there are no adjustable parameters whatsoever in 
Fig. |8j and that for ITO only the thickness was fltted. These results allow us to estimate 
the Casimir force expected in the two cases (Au and ITO) and compare the calculation with 



our measurements (see Eq. 12). 

The theoretical Casimir interaction is calculated with the Lifshitz equation using the 
dielectric properties of our surfaces. For Au, we have extrapolated the data of P^] with a 
Drude model {ujp = 9.0 eV and 1/r = 0.035 eV from [E]). For ITO, we used the model 
from for all frequencies. The computed force gradient is plotted as the black lines in 
Fig. [7j The agreement between data and theory is reasonable, although we do seem to 
obtain different powers for data and theory. At small separation, the experimental curves 
are bending upwards, which is a sign of surface roughness effects [13] ; the theoretical curves 
were calculated for perfectly smooth surfaces. Furthermore, the Au-Au data tend to give 
rise to a stronger force at large distance compared to theory, which is most likely caused 
by an artifact common to many AFM force measurements: the laser light is reflected from 
the planar sample into the photodetector giving rise to a background signal. In the Casimir 
force gradient method presented here, this artifact results, in flrst order, to an offset in the 
data; this explains the upwards trend of the data for large d. Although the precise distance 
dependence and strength of this artifact is unknown, we estimate from the force gradient 
data at large separation {d > 500 nm) that the associated systematic error is certainly 
< 2 N/m^. In the case of the Au-ITO measurements, such a background signal is a lot 
smaller because ITO does not reflect well the laser light (see Fig. [s] at w = 1.9 eV). To 
explain the mismatch, it is therefore more likely that the model we used for the dielectric 
properties [T3| is too metallic at low energy and that, consequently, the calculated Casimir 
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interaction is too strong especially at large d. 

So far, we have neglected the effects of surface roughness in our analysis. In Fig. |9j we 
show topology measurement of our surfaces. Figs. |9] a and b are AFM tapping-mode scans 
(10 by 10 fim) of the Au on polished sapphire and ITO on float glass samples, respectively. 
The gold sample has an RMS surface roughness of 0.8 nm, while the ITO coated plate has 
a roughness of 4 nm RMS. Fig. |9] c presents a height profile of the surface of the sphere 
bottom obtained with an optical profiler. Since the cantilever is mounted at a 15 degrees 
angle with respect to the planar sample surface (this is typical in AFM design), the top 
of this profile does not correspond to the area of closest approach in a force measurement. 
However, this height profile does give us the ability to estimate the surface roughness of the 
sphere, resulting in a value of 3.8 nm RMS. 

Since we used the same sphere in both sets of measurements, the surface roughness of the 
sphere can never cause the observed difference in Casimir force gradients between the Au- 
Au and Au-ITO cases. Furthermore, we recall that surface roughness tends to enhance the 
strength of the Casimir interaction [TB]. It is therefore impossible that the different surface 
roughnesses of the two planar samples is responsible for the difference reported in Fig. [7], 
because the ITO sample is considerably rougher than the gold coated sapphire substrate. 

When we discussed the experimental details of our experiment, we mentioned the inter- 
esting feature that we can measure both the Casimir force gradient and the hydrodynamic 



force acting on the sphere with the same lock- in amplifier at uj2- Fig. 10 shows the hydrody- 
namic force for both the Au-Au interaction (green triangles) and for the Au-ITO case (red 
squares). We have plotted the RMS force resulting from the 2.72 nm RMS oscillation of the 
plate at 002/ 2tt = 119 Hz. Both curves appear to change exponent at a separation of around 
200 nm. This bending is caused by the slip of the air flow across the surfaces, i. e. the fluid 
velocity at the gas-solid interface is nonzero. This phenomenon is treated in [16] and the 
expressions derived in there describe our data satisfactorily. Concerning the comparison of 
the two sets of hydrodynamic data, it is clear that the hydrodynamic forces are very similar 
in the Au-Au and Au-ITO experiments. Still, there exists a small difference between the 
two curves of roughly 2%. This difference cannot be caused by an error in the determination 
of the initial separation d^, because both data sets are parallel on the double logarithmic 
plot. We suppose that the cause may lie in the different surface roughnesses of the Au and 
ITO samples that lead to different amounts of fluid slip over the sample surfaces. 
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It is worthwhile to compare our method for the detection of hydrodynamic forces with 
recent measurements obtained with AFM's [71 [17]. In [7], the cantilever with sphere is driven 
at its free resonance and the amplitude and phase of the cantilever motion are used to extract 
the hydrodynamic force. In [17J, two methods were employed to measure the hydrodynamic 
interaction between a colloid sphere and a plate: measuring the static deflection of the 
cantilever during a fast approach of the planar sample and analyzing the thermal noise 
of the cantilever while slowly approaching the plate towards the sphere. In both papers, 
however, the separation between the two interacting surfaces was determined by bringing 
the sphere and plate into contact, a method that is prone to inaccuracies due to surfaces 
asperities (this is also reported in [17]). Since our method employs both a hydrodynamic 
force measurement and a precise calibration of the distance at the same time, we have 
developed a more reliable technique for hydrodynamic force measurements. 

CONCLUSIONS 

We have presented the experimental details of our Casimir force measurements between 
gold and ITO surfaces We have shown that the mechanical drift of our setup is less than 
0.1 nm per measurement run and that our electrostatic calibration is performed with 0.2% 
precision. Force gradient data obtained over approximately 72 hours reveal no drift in the 
signal at all, confirming the high stability of the setup. Furthermore, we have introduced 
and demonstrated a new method to determine the spring constant of our cantilever and 
the deflection sensitivity of the AFM readout. We also presented our measurements of the 
Casimir and hydrodynamic interactions between the gold and ITO surfaces, and provided a 
complete characterization of our samples in terms of their dielectric properties and surface 
roughness. 
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FIG. 1: (Color online) (a) Drawing of the experimental setup used to perform precise measurements 

of the Casimir force between a 100 ^m radius sphere and a plate. The aluminum block acts as 

a heat reservoir to keep the temperature of the setup constant. The instrument is based on a 

commercial AFM head that is, together with a custom-designed mechanical translator, mounted 

on the aluminum block, (b) Schematic representation of the working principle of the experimental 

technique. The Vb feedback circuit allows one to measure and compensate the residual voltage 

present between the sphere and the plate. The calibration lock-in amplifier is used to calibrate the 
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instrument and to find the initial separation between the two surfaces do. The measurement lock-in 




-2.0 -1.5 -1.0 -0.5 0,0 05 1.0 1.5 2.0 



tlo-<i^h("'^) 



FIG. 2: Mechanical stability of the experimental setup, (a) Mechanical drift in the initial separation 
do as a function of run number for all 580 Au-Au measurement runs. The error bars are determined 
by propagation of the error on a into the estimate of do by the fit with Eq. [8j The grey line 
represents a trend line that accounts for the slow thermal drift of the setup, (b) Histogram of the 
differences between the measured do values and the grey line of a. The line represents the best 
Gaussian fit, resulting in a 0.5 nm standard deviation. 
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FIG. 3: Stability of the electrostatic calibration, (a) All 580 obtained values for the force calibration 
constant k as a function of run number. The error bars are calculated by propagating the error 
on a. The grey line is a smooth trend line that accounts for slow variations, (b) Histogram of the 
relative deviations of the k data from the trend line in c. The Gaussian fit has a standard deviation 
of 0.2%. 
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FIG. 4: (a) The data points represent measurements of the total force gradient as a function of 

separation between the sphere and plate surfaces for 150 measurement runs out of a total of 580 

runs. The line shows the electrostatic force gradient associated to the simultaneous calibration 

procedure, (b) Plot of all 580 force gradient measurements obtained for d ~ 95 nm as a function 

of time. The grey line represents the average force gradient, (c) Histogram of all the relative 

deviations between the single force gradient measurements around 95 nm and the average force 

gradient. The Gaussian fit has a standard deviation of 3.5%. 
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FIG. 5: Plot of the electrostatic force gradient difference between a measurement run performed 
with a strong electrostatic interaction and a run performed with a weak electrostatic force. Data 
are plotted as a function of separation. The line corresponds to the electrostatic force gradient 
obtained from the calibration signal. See text for details. 
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FIG. 6: Plot of Vac \/ S2uji 1 5'4wi as a function of piezoelectric transducer extension. The line 
represents the best fit of the data with Eq. 19 (reduced = 0.25). The slope of the fit can be 
used to obtain the cantilever spring constant and the deflection sensitivity 7. 
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FIG. 7: (Color online) (a) Casimir force gradient as a function of separation for the Au-Au (green 

triangles) and Au-ITO (red squares) interactions for randomly chosen subsets of the data (150 out of 

580 for both cases) plotted on a double logarithmic scale, with the common electrostatic background 

subtracted from the data. The black lines correspond to the calculated Casimir interactions. 

(b) Histogram of all 580 force measurements for both the Au-Au and Au-ITO measurements at 

d = 120 nm. The difference in Casimir force gradient is ~ 50% between the gold and ITO 

measurements, (c) Same as b, but for d = 80 nm. At this separation, the difference in the force 
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gradient is ~ 40%. 
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FIG. 8: (Color online) Measured reflection (R) and transmission (T) spectra as a function of 
photon energy for the Au on sapphire (green data) and ITO on float glass (red data) samples. The 
continuous black lines are calculations of the reflection and transmission spectra expected for our 
samples (no fit parameters except ITO layer thickness), using handbook data for gold [12] and a 
model from [13] for ITO. The transmission spectra for gold are zoomed in because the maximum 
transmission (around 2.5 eV) is only 1.4%. The calculation of the transmission spectrum of the 
ITO sample is quite sensitive to the choice of dielectric properties of the float glass for photon 
energies above 4 eV. The black lines describe the measured data reasonably enough to allow for 




FIG. 9: (Color online) (a,b) AFM topography scans (10 by 10 fim) of the surfaces of the Au on 
sapphire and ITO on float glass samples, respectively. The surface roughness for the gold sample 
is 0.8 nm RMS, while the ITO plate has a surface roughness of 4 nm RMS. (c) Optical profiler 
scan of the bottom of the gold coated polystyrene sphere that is attached to the cantilever for our 
force measurements. The surface roughness of the sphere is 3.8 nm RMS. 
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FIG. 10: (Color online) Hydrodynamic force (RMS) acting on the sphere as a function of separation 
for a subset of the data (150 out of 580 for both cases), caused by the oscillations of the plate surfaces 
at 119 Hz. The green triangles represent the force in the case of two gold-coated surfaces, while 
the red squares correspond to the Au-ITO interaction. 
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